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Abstract 

The extend to which geometrical effects contribute to the production and sup- 
pression of the J/V' and qq minijet pairs in general is investigated for high energy 
heavy ion collisions at SPS, RHIC and LHC energies. For the energy range under 
investigation, the geometrical effects referred to are shadowing and anti-shadowing, 
respectively. Due to those effects, the parton distributions in nuclei deviate from 
the naive extrapolation from the free nucleon result; /a 7^ Af^. The strength of 
the shadowing/anti-shadowing effect increases with the mass number. Therefore it 
is interesting to see the difference between cross sections for e.g. S + U ys. Pb + Pb 
at SPS. The recent NA50 results for the survival probability of produced J/'i/^'s 
has attracted great attention and are often interpreted as a signature of a quark 
gluon plasma. This publication will present a fresh look on hard QCD effects for 
the charmonium production level. It is shown that the apparent suppression of 
J/ip^s must also be linked to the production process. Due to the uncertainty in 
the shadowing of gluons the suppression of charmonium states might not give reli- 
able information on a created plasma phase at the collider energies soon available. 
The consequences of shadowing effects for the xp distribution of J at ^/s = 20 
GeV, -y/i = 200 GeV and ^/s = 6 TeV are calculated for some relevant combina- 
tions of nuclei, as well as the pT distribution of minijets at midrapidity for A'^^ = 4 
in the final state. 
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1. Introduction 

Since the advent of QCD in the 70's great emphasis was laid on the 
existence of a phase transition of, yet unknown, order, being typical for non- 
abehan gauge field theories. From lattice calculations it was emphasized 
that, at zero chemical potential, a phase transition should show up at some 
temperature Tc ~ 150 — 200 MeV when explicitly taking quarks into account. 
The value for Tc is slightly higher for a pure gauge theory. Also, at non-zero 
chemical potential, as suggested in the MIT bag model, one should access a 
phase transition due to the increasing outward pressure of the partons inside 
the bag finally leading to a deconfined phase. Due to the difficulties emerging 
when considering dynamical fermions the work on non-zero chemical poten- 
tial has not yet reached the same level of success as that for /x = in lattice 
QCD. 

Now, in actual high energy heavy ion collisions the following scenario can 
occur. Two streams of initially cold nuclear matter collide and may result 
in a plasma phase, which is created within the transverse dimension of ap- 
proximately the size of the overlapping nuclei. The plasma cools down to 
form hadronic degrees of freedom in the subsequent expansion. If one has 
this phase transition in mind one also has to confront the question of its 
experimenatal detection. Typical signatures under discussion are leptonic 
(dilepton [0] and photon production due to the interactions among the 
quasi free partons via the different QCD processes qq — > 7(7, gq 7g, ...) 
and hadronic ones, such as the suppression of J/ip^s. Now, the QCD reac- 
tions in the plasma are not the only source for leptons. One expects a large 
background coming from the decay of vr" and rj mesons. It therefore is nec- 
essary to carefully handle this background by experimental methods such as 
invariant mass analysis. 

It is also obvious that the signatures have to give clear and powerful in- 
formation on the plasma phase. Escpecially when looking at the hadronic 
signatures this does not have to be the emphasized in where it was 

shown that gluon depletion due to DGLAP splitting in the colliding nuclei 
can lead to the same results as the current experiments at NA50 show, 
which in turn implies that those experimental results propably have lost their 
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meaning as a plasma signature at SPS. 

2. Production and suppression of the J/ip 

The J/ip is a cc bound state interacting via two forces in a confined sur- 
rounding: a hnear confining potential and a color-Coulomb interaction. In 
the plasma phase the linear potential is absent due to the high temperature 
leading to deconfinement. Every color-charge is Debye-screened by a cloud of 
surrounding quark-antiquark pairs which weaken the binding force between 
the cc pair, thus reducing the color charge seen by the other (anti) quark. 
Since the density of the screening pairs rises strongly with increasing tem- 
perature, the binding force gets weaker and weaker when the temperature 
rises above Tc. As a result, the charm quark and antiquark drift away from 
each other, so that finally no bound state formation is possible in a plasma 
phase of high enough temperature [§] . 

However, the plasma phase is not the only source of suppression 0. One also 
has to take into account final state interactions for this hadronic degree of 
freedom that are absent for leptonic signatures. Because the J /if) is a, verly 
weakly bound state, the interaction with nucleons and secondaries, that are 
always present in a heavy ion collision, in addition significantly lowers the 
survival propability for a J It is obvious that such effects should increase 
with increasing mass number. One also expects the phase transition to hap- 
pen for the heavier nuclei. Therefore one has two effects both increasing with 
the number of nucleons involved. This in turn implies that the experiments 
have to be done with very high precision to disentangle those effects. 
At this point another source of suppression comes into play that also in- 
creases with the mass number and therefore has to be accounted for: nuclear 
shadowing. This effect already enters on the production level of the charmo- 
nium bound state. The former two effects, namely suppression by melting in 
the plasma phase and comover activity, enter only at a level when the J/ip 
already exists at later proper times r. Now the nuclear shadowing effect ap- 
pears when the charmonium is produced via the various processes depicted 
in figure 0. 

The total hidden charm cross section in pN collisions below the open charm 
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Figure 1: The various LO processes leading to the direct production of a cc 
pair. 

threshhold is given by 

CTccis) = ds dXAdXBfi{xA)fj{xB)^{s)S{s-XAXBs) (1) 

J4m,2 J 

Here, fi and fj denote the parton densities and a is the cross section on the 
parton level, i.e. qq cc, gg — >■ cc. The cc pair subsequently will turn into 
a color singlet by interaction with the color field, induced in the scattering, 
the so-called " color-evaporation" mechanism. In the J/ ip production in a 
proton nucleon reaction was parametrized as 

= fj,^^%''' (2) 

with f^i^ = 0.025 from comparison with data Here the production of 
the J/ip is described as proceeding via the NLO production of a ccg state 
and a subsequent evaporation of the gluon. (For a more detailed model 
also including the non-relativistic quarkonium model in the quarkonium- and 
bottonium- nucleon cross section see 0). 

It is obvious that any changes in the parton densities will result in changes 
of the cc production cross section. Because we know since the EMC mea- 
surements ||I0| that //^ 7^ Afl^ this demands for some further investigation. 
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Here we will investigate the influence of the nuclear gluon and quark distri- 
butions on the J/iIj production cross section by using a modified version of 
a parametrization based on a (impact parameter averaged) data fit given in 
in]] . We will show the influence on the differential cross section da^^ /dxp 
for gg fusion and qq annihilation given by 

X g\xA,QV{xB,Q') (3) 



X 



[q^{xA,Q')q''ixB,Q')+q^q\ (4) 
and on the minijet cross section 

da A da^^^'^^ 

2'KXAg {xA,PT)xBg^{xB,pD (5) 



PTdpTdyidy2 dt 

at midrapidity y = yi = 1/2 = 0. We choose rUc = 1.5 GeV and mn = 
1.85 GeV. The momentum fractions are given as xa,b = 1/2 [ixj? + (x^^ + 
4Q^/s)^/^] for the xp distribution and x = 2pt/ \/s for the minijets at midra- 
pidity. We take all cross sections in leading order and do not include any K 
factor for higher order contributions since we are mainly interested in relative 
effects. For the parton distributions we choose the CTEQ4L parametriza- 
tion. 

The reason for our investigation is the following: the recent NA50 data show 
a deviation from the tendency expected from earlier experiments when the 
mass number of the involved nuclei is increased. Now, in it was shown that 
due to multiple scatterings between partons the uncertainty in the survival 
propability gets so large that one cannot distinguish whether the data found 
by NA50 is due to gluon splitting in the production phase or due to plasma 
absorption as claimed by several authors. Obviously, the originally good idea 
of J/ ip suppression as a good tool for plasma investigation seems to has lost 
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its predictive power at the available energies. It is therefore interesting to 
see what one can expect at future colliders. 

In the next part we will give some details of the parton densities in nuclei 
in the energy regimes of SPS, RHIC, and LHC. 



3. Nuclear shadowing and the connection to the J/il^ 

The history of the modification of nuclear structure functions, as com- 
pared to the free nucleon ones, is founded on the findings of the EMC group 
that lead to the so-called EMC effect |Ty] (even though one should say that 
shadowing effects in principle have been known since the 70 's This 
effect shows that fA 7^ ^/at, which implies that the parton density in the nu- 
cleus is not simply given by the nucleon number times the respective parton 
density in the nucleon. Depending on the frame (lab- or infinite momentum 
frame) one derives completely different interpretations for the nuclear struc- 
ture functions and for the deviations from the naive pp-extrapolations. For 
typical values of the momentum transfer in a pp reaction of pt = 1 — 6 GeV, 
where perturbation theory should be applicable, one is in the so-called anti- 
shadowing region for SPS and in the shadowing region for RHIC and LHC 
at midrapidity. In the following we will shortly review the interpretation 
of shadowing in the two relevant frames and will start with the lab frame 
description which is the natural frame for typical deep inelastic scattering 
measurements off nuclei (at least from the experimental setup point of view). 



A. Lab frame description 
In the lab frame the expression shadowing immediately seems to imply 
a geometrical effect. When one speaks of something lying in the shadow of 
another thing one means that the second body is not visible since the first 
body is placed nearer, e.g. to some source of light. A similar reasoning can be 
applied in the case when a lepton is scattered off a nucleus consisting of many 
nucleons. The exchanged virtual photon does not (in the relevant x range) 
interact individually with each nucleon but coherently with all nucleons or 
at least with a major part of the nucleons inside the nucleus; some nucleons 
are therefore lying in the shadow of other (surface) nucleons. As we will 
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Figure 2: The two physical processes arising from the two possible time 
orderings. 

see later, this reasoning is linked to the momentum fraction of the struck 
parton inside the nucleon. Because the momentum fraction is bound from 
above this interpretation is limited to the explanation of shadowing and is not 
applicable for the reasoning of anti-shadowing, the EMC effect or the Fermi- 
motion effect. Unfortunately there is yet no single theory to understand the 
whole range of the momentum fraction from < x^j. < 1. For an excellent 
review of different models and interpretations see [|T4| . 



For a deep inelastic scattering process there exist two possible time orderings 
for the interaction of a virtual photon with a nucleon or with a nucleus: either 
the photon hits a quark inside the target (the so-called hand-bag graph) or 
the photon creates a qq pair which then strongly interacts with the target. 
Those two possible processes are depicted in figure ^. 

As can be seen from the ratio of the amplitudes of the two processes one 
realizes that the diagram on the right hand side only contributes at small 
enough a; (x <C 0.1). At low the interaction of the virtual photon with 
the nucleons inside the nucleus happens via the low mass vector mesons p, uj 
and as described in the vector meson dominance model (VMD) with the 
typical spectral ansatz for the description of the fluctuation spectrum [|I^]. 
The reduction in the quark density, manifesting itself in the shadowing ratio 
= /A ■ , can then be understood in terms of a multiple scattering 
series where the fluctuation interacts with more than one nucleon over a 
coherence length of Ic ~ l/{2mx). At higher the partonic degrees of 
freedom are probed; nevertheless, shadowing is due to long distance effects 
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and therefore always incorporates a strong non-perturbative component, even 
at large . Also, the qq continuum has to be taken in addition to the mesons 
giving rise to the generalized VMD model. 

The interaction of the virtual photon with a nucleon can essentially be split 
up into two parts: the virtual photon with its quark- ant iquark fluctuation 
and the interaction of the fluctuation with the parton which happens via 
gluon exchange: 

a(7*iV) = J d\\^{z,r)\'aggN{r) (6) 

where the Sudakov variable z gives the momentum fraction carried by the 
quark (or antiquark). 

The cross section for the interaction of the fluctuation with the nucleon 
can be described in the DLA as 

a,,N = ^r'a.{Q")x'g{x',Q") (7) 

where x' = Mgg/(2mz/), r is the transverse separation of the pair and Q'^ = 
Due to 

mz,r)\^^ (8) 

pairs with small transverse separation are favored. As can be seen from (|^) 
this in turn implies a small cross section. This smallness is compensated by 
the strong scaling violation of the gluon distribution in the small x region 
as r (Q'^) decreases (increases). In the Glauber eikonal approximation the 
interaction with the nucleus is expressed in terms of the nuclear thickness 
funtion TA^b) as 

a^^A = I d% (l - e-'^---^-^(^)/2) (9) 

When there is a longitudinal momentum transfer appropriate to the produc- 
tion of the hadronic fluctuation h a phase shift behind the target results and 
the incident wave exp{ik2z) is changed to —T {b)ex]i{ik'2 z) with kj 7^ fc^ and 
the nucleus profile function r(6). The phase shift Ak^ = kj — k^ in turn gives 
rise to the coherence length ~ l/^^^- When now ^ 2Ra (i.e. for small 
momentum transfers Akz), the hadronic fluctuation interacts coherently with 
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Figure 3: Graphic representation of the lab frame interpretation of shadow- 
ing. 



all nucleons inside the nucleus, Glauber theory is valid and a reduction in 
the cross section results (for further details we refer to |]I4|, ^ |l^, 0). For 
an illustration of the effect see figure |^. 

B. Infinite momentum frame description 
In the infinite momentum frame a completely different mechanism is em- 
ployed. The key idea here is the fusion of partons giving rise to a process 
that competes with parton splitting expressed in the DGLAP equations. This 
idea was first formulated in and later proven in [|l8l . In the following we 



will give the main ideas and conclusions of the parton fusion model [19 



As is known, in the infinite momentum frame the Bjorken variable xbj is in- 
terpreted as the momentum fraction of a parton with respect to the mother 
nucleon. When now, inside a nucleus, the longitudinal wavelength of a par- 
ton exceeds the Lorentz-contracted size of a nucleon or the inter-nucleon dis- 
tance 2i?7v, then partons originating from different nucleons can "leak out" 
and fuse. This effect can be estimated from l/{xP) ~ 2RnMn / P to show 
up at X values smaller than x ~ 0.1. As a result of the parton-parton fusion 
partons are "taken away" at smaller values of x and "shoveled" to larger 
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values of x where anti-shadowing appears to guarantee momentum conserva- 
tion. As a result of the parton fusion the x-range for the measured structure 
function is expanded to values x > 1. Hereby, an alternative description 
of Fermi-motion is achieved. In the lab frame interpretation the saturation 
of shadowing was interpreted in terms of a coherence length larger than the 
nucleus. Here, the saturation towards smaller x values is interpreted in terms 
of the longitudinal parton wave length exceeding the size of the nucleus. In 
that sense the infinite momentum frame interpretation of shadowing is for- 
mulated in terms of variables that are inherent to the nucleus and there is 
no need for a scattered lepton or a collision. 

In addition to the longitudinal shadowing one expects an additional shadow- 
ing effect from the transverse fusion of partons: for sufficiantly small values 
of X and/or Q"^ the total transverse occupied area of the partons becomes 
larger than the transverse area of the nucleon. This happens (e.g. for gluons) 
when xg{x) > Q^B? where the transverse size of a parton is l/Q^ and R 
is the nucleon radius. The depletion in the gluon and sea-quark densities 
arising from that process are expected at values x < 0.01. 

4. The used parametrization 
In ig a fit to the E772 [0, NMC |21[ and SLAG data was given 
as a parametrization for the ratio Rf2 = F-i 1^ ' ^2 (see figure ^). Now 
this parametrized ratio cannot simply be multiplied with all the individual 
parton distributions entering the formulas. One has to make a distinction 
between the valence and the sea quarks and also needs a different ratio for 
the gluons. Our results for da'^^/dxp are based on the shape of the ratio 
given in ||ll| at the initial momentum transfer Qq = 2 GeV. Up to now, the 
production processes were often calculated by using the measured shadowing 
ratio Rp^. From the lab frame interpretation we know that the cross section 
for the interaction of a gluon pair is larger than the one for the interaction 
of the quark-antiquark pair (Cgg^v = S/^Cgliv )• The same tendency can be 
found in the parton fusion model. In calculations in the parton fusion 
model for -^^^Sn showed an impact parameter averaged gluon shadowing that 
is as twice as strong {Rq ~ 0.34) as the sea quark shadowing at x = 10~^ 
and = 5 GeV^ already for this light nucleus. To account for the much 
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Figure 4: Fit to the data for various nuclei at Qo = 2 GeV as given by Eskola. 
stronger gluon shadowing we therefore modified the parametrization given in 
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In the lab frame, the relevant range for the coherence length to produce 
the shadowing effect is Iq = rj^^ ~ 1-8 fm < Ic = l/(2mx) < 2Ra. For 
Ic ^ '^Rai corresponding to x <^ 0.1 fm/ 1.1 fmA^^^, the shadowing of 



gluons at some initial scale at fixed impact parameter behaves as p5| 

Aeff _ 2 - 2(exp - R/2) 



A R 



(10) 



where R = T{h) ■ <Jeff- For the interaction of the qq pair one finds creff,qq ~ 
14 mb. which approximately corresponds to the pN cross section. We here as- 
sume that the perturbative factor 9/4 is valid also for the non-perturbative 
regime and therefore choose (Teff,gg ~ 30 mb. At 6 = and for Pb one 
therefore has a maximum amount of shadowing of A^ff/A 0.39 which is 
approximately 15% smaller than the 6-averaged result. Because the two dif- 
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Q' - 10 GeV' 

1.2 - Q' = 4GeV' 




10"* 2 5 10"' 2 5 10"^ 2 5 10"' 2 5 10° 

X 



Figure 6: Gluon shadowing ratio evolved to = 10 GeV^ with DGLAP 
without fusion terms. Due to the narrow range 4m^ < < Amj^ in the 
interpretation we use the ratio at some fixed intermediate scale = 10 
GeV^. 
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ferent scenarios (lab- or infinite momentum frame) give sucli different results, 
Rg ~ 0.39 |2S[ vs. Rg -C 0.3 for heavy nuclei at small x, we decided to 
choose some intermediate value as a starting point for the DGLAP evolution. 
We therefore employ the curves shown in figure | to account for the large 
difference in the quark- and gluon shadowing ratios. Due to the large uncer- 
tainty of the initial Rg we choose the same ratio for Au and Ph at = 10 
GeV^ as shown in figure ^. Also, we again want to emphasize that the com- 
monly used shadowing ratios, only account for impact parameter averaged 
measurements in DIS reactions. Therefore our results should be seen as for 
central events only because the production mechanism in very peripheral 
collisions should produce significantly smaller rates with significantly smaller 
influences from shadowing effects p6| . 

For the minijet cross section we used a Q^-dependent parametrization given 
in 1 16] to account for the larger region. 

4. Results 

We will first present the results for the minijet cross section including only 
processes i,j — > k, I with i,j = g and k,l = qq with four flavors in the final 
channel (due to the dominance of the gg fusion process annihilation processes 
are neglected at RHIC and LHC at midrapidity) . For the xp distribution 
we used our modified version of the parametrization in |TT| but for the mini- 



jets we used an impact parameter dependent parametrization with b =0 [|16 



shown in figure |^. This parametrization is applicable here since we are in the 
pure shadowing region where the generalized VMD approach used to derive 
it is applicable (even though one should say that the Glauber ansatz should 
only be valid up to values x ~ 10~^ as restricted by the eikonal approxima- 
tion). 

The regions of the momentum fractions corresponding to the momentum 
range 1 GeV < pr < 6 GeV for RHIC (v^ = 200 GeV) and LHC (v^ = 6 
TeV) are represented in figure |^ as shaded areas. The results for the cross 
sections for RHIC and LHC are given in figure |^. In that calculation all 
quark antiquark pairs {k, I = qq) up to the bottom threshhold were taken 
into account, i.e. A'^^ = 4 in the final state. One clearly sees the deviation 
having its origin in the shadowing of the nuclear parton distribution. As ex- 
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Figure 7: Gluon ratios corresponding to the various energy and transverse 
momentum regimes. 
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pected, the shadowing effect decreases as pr increases due to the momentum 
fraction x = 2pt/ \fs. 



Next we will present the results for the da^^/dxp cross sections. We 
first calculated the proton-proton cross sections to show the dominance of 
the gg fusion process over the qq annihilation process at small xp (see figure 
^): The results for the cross sections for S + U and Pb + Pb at SPS at 
y/s = 20 GeV, Au + Au at RHIC at = 200 GeV, and Pb + Pb at LHC 
at = 6 TeV are presented in figure 0. In this case we were restricted 
to our modified version of the parametrization of |jll| due to the integration 
reaching up to momentum fractions x > 0.1, not allowing us to use the same 
parametrization as for the minijet production. At SPS energies one clearly 
sees the different regions of the parametrization entering the cross section. 
At small xp one has the enhancement due to the antishadowing which is 
followed by the depletion due to the EMC region at larger x and finally 
one can identify the Fermi motion effect as — > 1. The effects are clearly 
stronger for Pb + Pb than for S + U (compare figure To get an impression 
of the relative strength of the nuclear modifications in the respective nuclei we 
calculated the ratios of the shadowed to unshadowed cross sections in figure 
|TT|. The difference between Pb + Pb and S + U at SPS energies in principle 
is only small; in the relevant region of small xp, where the cross section has 
not dropped yet too much, the charmonium production in Pb + Pb is slightly 
larger than in S + U (~ 8%). For RHIC one is in the shadowing region. 
The suppression strongly varies over the xp range between ^ 0.6 — 0.35. At 
LHC an even stronger suppression is found due to the smaller momentum 
fractions entering the shadowing ratios. Here the suppression is ~ 0.3 — 0.5. 
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;ure 8: Minijet cross sections for creation of qq pairs for RHIC and LHC. 



17 




18 




19 




Figure 11: -^^i^/^^^"^^ for SPS (P6 + Pb,S + U), RHIC (Aii + AC/) 
and LRG (Pb + Pb). 
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5. Conclusions 

From the results shown above one now can draw the following conclusions 
for the consequences of the shadowing effects for charmonium production and 
suppression at SPS, RHIC and LHC. 

First, one can conclude from figures ^ and that an enhancement of char- 
monium states produced near midrapidity due to antishadowing at ^/s = 20 
GeV is predicted (small xp). For larger xp, a clear suppression of the charm 
cross section to ~ 70 — 80% of the unshadowed result (figure |1^ and again 
a rise at the largest xp values is predicted (the latter one due to the Fermi 
motion effect). 

For RHIC energies of y/s = 200 GeV the situation changes; for minijets with 
1 GeV < pt < 6 GeV at midrapidity (or at small xp, respectively) one is 
completely in the shadowing region. Here, the shadowed result are reduced 
by ~ 45%. At LHC the situation is even more dramatic: the ratio of the 
shadowed cross section to the unshadowed cross section at p^- = 1 GeV is 
0.22 which amounts to a suppression of a factor ^ 4.6. 
Similar effects are observable for da^^^'^^'^ /dxp: at small xp ~ 0.05 for 
RHIC the cross section is reduced by a factor ^(jshad. j ^^unshad. ^ o.58, and 
gets suppressed even more towards larger xp down to values ~ 0.35. At 
LHC one finds a less strong variation over the xp range with a mean value 
of ~ 0.35. In these results one problem is unveiled: the difference between 
the (not yet exactly known) gluon ratio Rq and the quark ratio Rp^ that, 
according to the calulations in |^ increases with increasing mass number. 
If, as it was recently done at CERN-SPS, the future experiments at RHIC 
and LHC compare different combinations of nuclei and derive results similar 
to the NA50 data one has to ask oneself whether one has detected the plasma 
or whether the detection is that the gluon ratio in not simply given by Rp2, 
even at small x. To give clear predictions it is mandatory to control the 
value of Rg at the typical semihard scale Qsh GeV with high precision. 
Therefore charmonium and bottonium suppression effects can also be due to 
purely geometrical effects, i.e. shadowing. 
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